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A comprehensive study of doped RVB states is performed. It reveals a fundamental connection 
between superconductivity and quantum spin fluctuations in underdoped cuprates : Cooper pair 
hopping strongly reduces the local magnetization mo. This effect pertains to recent muon spin 
rotation measurements in which mo varies weakly with hole doping in the poorly conducting regime, 
but drops precipitously above the onset of superconductivity. Gutzwiller mean field Approximation 
(GA) is found to agree with numerical Monte Carlo calculation. GA shows for example that for a 
bond amplitude u{r) = e~ r/ % spin spin correlations decay exponentially with a correlation length 
oc e 3 ^ 4 //2 . Expectation value of the Heisenberg model is found to be correlated with average loop 
density. 



I. INTRODUCTION 



When holes are introduced into the copper oxide planes of high T c Cuprates, spin and charge correlations change 
dramatically. The local magnetization mo, measured by /iSR B on e.g. L^-rrSr^CuC^, reveals a qualitative difference 
between the insulating and superconducting phases: mo is rather insensitive to doping in the poorly conducting regime 
< x < 0.06, but drops precipitously above the onset of superconductivity at x > 0.06, becoming undetectable at 
optimal doping x ~ 0.15. Theoretically, holes can cause dilution and frustration in the Heisenberg antiferromagnet, 
which create spin textures: either random ("spin glass") or with ordering wavevector away from (it, it) (sometimes 
called "stripes") ||. However, the apparent reduction of local magnetization by the onset of superonductivity, is a 
novel and poorly understood effect. Theory must go beyond purely magnetic models, and involve the superconducting 
degrees of freedom. 

We find that this problem is amenable to a variational approach, using hole-doped Resonating Valence Bonds (RVB) 
states. The RVB states were originally suggested by Anderson to describe the spin and charge correlations in the high 
T c Cuprates §. They are excellent trial wave functions for the doped Mott insulators, with large Hubbard repulsion 
U since: 

(i) Configurations with doubly occupied sites are excluded. 

(ii) Marshall's sign criterion for the magnetic energyj3jis satisfied, and Heisenberg ground state energy and antifcr- 
romagnetism at zero doping is accurately recovered |7H. 

The hole-doped RVB state is a new class of variational states, in which spin and charge correlations are parameterized 
independently, without explicit spin nor gauge symmetry breaking. Such parameterization allows states with magnetic 
and independently d or s-wave superconducting (off-diagonal) order or disorder, thus permit an unbiased determination 
of ground state spin and charge correlations appropriate for the Cuprates. These are important advantages over 
commonly used Spin Density Wave, Hartree-Fock and BCS wavefunctions. 

A comprehensive study of the state is performed using Monte Carlo and mean field calculations. Phenomenological 
low energy effective Hamiltonian is proposed, with two major components: Heisenberg interaction for spins, and single 
or Cooper pair hopping kinetic energy for fermion holes. 

Regarding this model our key results are: 

(i) For the magnetic energy alone, the local magnetization mo is weakly dependent on doping concentration. This 
holds independently of inter-hole correlations for either randomly localized or extended states. 

(ii) In contrast to (i), mo is strongly reduced by the kinetic energy of Cooper pair hopping, which correlates the 
reduction of mo with the rise of superconducting stiffness, and hence j^] the transition temperature T c . 

These results agree with the experimentally reported correlation between mo and T c jD]. This relation appears to 
be weakly dependent on the precise hole density. 
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We also find that RVB states have the following properties: 

(i) The magnetic energy is correlated with the average loop density: 

r = L 2 to§/ (average radius of gyration of a loop) 2 , where L is the linear size off the lattice. 

(ii) The Gutzwillcr mean field Approximation (GA) for magnetic correlations is in good agreement with the Monte 
Carlo results. 

(iii) Long range magnetic correlations in RVB states are extremely sensitive to changes in the singlet bond amplitude 
u. For example with u(r) = exp (—»"/£) the spin-spin correlation function decays exponentially with correlation length 
£ ex oc exp ((1 — x)^^ 2 ), where x is the hole concentration. 

The paper is organized as followed: Sec.|| introduces the hole-doped RVB state, and discusses the numerical proce- 



dure. Sec.[llj defines our variational parameters. Sec. IV deals with the antiferromagnetic and superconducting order 
parameters. Sec|v| deals with the com ponents of the effective Hamiltonian. Sec. VI correlates between superconducting 



T c and local magnetization. Sec. |VIl| is a summary and discussion. 

The paper has 3 appendices. App. ^ reduces the hole part of the doped RVB to a numerically convenient format. 
App. [b| derives expressions for expectation values. Particularly, an alternative procedure to calculate magnetic 
correlation is derived and used to check the computer program. In App.^, the GA is performed analytically. 

II. THE HOLE DOPED RESONATING VALENCE BOND STATES 

A Valence bond (VB) state is 

i«>= n (oM-&M)i°>> 

(i,j)£a 

where a is a pair covering of the lattice, aj, &t are Schwinger bosons, and i — 1, . . . L 2 is a site index on a square 
lattice. 

RVB states are superposition of VB states. We restrict the disscussion to 

i*n)=e n - 6 W)i°>- ( 2 ) 

a (i,j)ea 

where u(r^ ) > is a variational singlet bond amplitude, which connects sites of different sublattices A and B only. 
This ensures Marshall's sign ||. 

The hole doped RVB state is defined by : 



|*[u,«;ar]) =V G {x)\i>[u,v\) 



\ip[u, v}} = exp 



E + Uii(ot6}-6t ]) 



10) (3) 



where fj are spinless hole fermions, Uij > 0, and v{vij) is an independent hole bond parameter. The Gutzwiller 
projector Vg{x) imposes two constraints. A constraint of no double occupancy: 

< + r4+n} = l Vi, (4) 

and a global constraint on the total number of holes: 

53 n} =xL 2 = N h . (5) 

i 

Due to Va{x), W can be written as a sum over bond configurations of singlets and hole pairs which cover the lattice 
as depicted in Fig. |l|. 

An overlap of two VB states, (cx\/3), is expressed in terms of a directed loop covering of the lattice (DLC) ||7|, ^0| , p^[ , 
hence : 

(*[u]|*M) = £q a (6) 

A 

where A is a DLC, 
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"a = n h n ^ ) ( ? ) 



and A is a directed loop. 

With the results of App.[A| the norm of the doped RVB state is 



($[u,v;x]mu,v;x})= ^ W( 7 ,A( 7 )), (8) 

7,A(7) 



where 7 is a distinct configuration of Nh holes sites: 

7 = {(it € A, j fe € S)} fe ^ 1 : Vfc < i fc < i fc >, j fe < jV, 
A 7 is a DLC which coveres the lattice but the hole sites, 



^(7,A 7 ) = { det 2 V(7)^ *>0 (g) 
and V is an A^/2 x matrix with 

V{l) kl = v ikjl . (10) 
Expectation value of an operator O is expressed as a weighted sum 

<°> = 7*T*\ E ^(7. A 7)0(7, A 7 ) = O (11) 

where 0( 7 , A 7 ) is defined by Eqs. @ an d (FT[) . 

We use standard Metropolis algorithm |jl 1|| for the evaluation of sum ([ill ) ■ The basic Monte Carlo step for updating 
the DLCs is the one used by Ref. : Choose at random a site and one of its next-nearest neighbors and exchange, 
with transition probability that satisfy detailed balance, the bonds connecting each of them, either to the next site 
(forward-bond), or the previous site in their loops. In Ref. [ fl2| we show, that for u r > Vr, these steps are ergodic, 
that is, any DLC can be reached from any other by a sequence of Monte Carlo steps. 

For the fermion holes our update scheme is a simple generalization of the "inverse-update" algorithm of Ceperley, 
Chester and Kalos (T^| . According to Eq. (|l0|) , changing a position of an A (B) sublattice hole amounts to changing 
one row (column) in the matrix V. In the our calculation boundary conditions are periodic. 

For dimer doped RVB state, where 



J 1 otherwise v ; 



we obtained exact results using transfer matrix technique. For a 4 x 40 undoped lattice fl4|| the magnetic energy 
is E mag = -0.320744 J /bond. The Monte Carlo result is -0.3210 ± 0.0002. (SoS r ) is exponentially decaying with 
correlation length of ^(jj mer = 0.724, the Monte Carlo result is 0.738. Exact and Monte Carlo results for the doped 
2 x 64 ladder jl2) appear in Fig.(||). Our program successfully reproduced existing data for RVB states f^j8|,[l^|. Other 
tests of the program appear below. 

We also use the Gutzwiller Approximation (GA) to evaluate expectation values in the doped RVB state. The GA 
is discussed in App.|c]. 

III. THE VARIATIONAL PARAMETERS 

In the undoped, x = 0, case we treat three classes of singlet bond amplitude u. 

Mr) = y P ( 13 ) 
u ex (r) = ii sr (r)^i exp (14) 
u g (r) = u sr (r) exp (-Qr 2 ) (15) 
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with u(l) = 1 and 

u sr (r) = ai exp + a 2> ( 16 ) 

where for u ex (u g ) £^ = 1.7 (2) and a 2 = 0.05 (0.018). u sr determines the short range decay of u ex and u „ p2| . We 
also use it — umf- umf is derived from the Schwinger-boson mean field theory of the Heisenberg model ]15|,HI- For 
x > we use u p , Eq. (^3|), and u e:E , Eq. (|l4|). 

For the function v the following cases of inter-hole correlations are treated: 



,( r u) 







„,met 1 / f2v^ met — ikr 

w «(r) = c a(^) (5 r : f ) , a = s,d 



sign(k) k G £ 
k g £ 



(17) 



where |w^ le | = 1, fj are nearest neighbor vectors on the square lattice, c s = 1 and ca = fj x — fc. 

v] ns puts the Nh holes on random sites. This state describes an insulator with disordered localized charges. 
v me t describes weakly interacting holes in a "metallic" state: 

n /ki°) =pg(*)«p fE< et /k/-k +( ^))i°) =^(a ; )cx P (E^rv/znio) (is) 

kes \ k / \ ij / 

where the product is over Nh states, 



(19) 



and v% et = X; k < et e _ik(r ^ rj) . Here we check S which is centered at k mi „ = (±f ,±f). See Fig. |. Results for £ 
centered at k m i„ = (0, ±7r) , (±tt, 0) are no£ qualitatively different ]l2| ). v met obey 

«k+(7r,7r) = -«k, (20) 

hence v^ et only connects i £ A to j S -B. Correlations in a state with v = u m et were previously computed by Bonesteel 
and Wilkins ||. 

u s and Vd describe tightly bound hole pairs in relative s and d-wave symmetry respectively. 



IV. ORDER PARAMETERS 



A. Local magnetic moment and long range magnetic correlations 

The local magnetization is 

mg = -^^(S,S J -)e- i ^)^- r ^ (21) 

ij 

where e.g. S+ = S x + iS y = atb. With respect to Eq. (|ll|), (SiSj) = S{r io ) is calculated using Q 

co/ A^ I if i j are on the same loop in , nn s 
S ^(T' A 7)i = ( otherwise (22) 

where the sign is + if i and j are on the same sublattice. To check our program we also used an alternative procedure 
to calculate magnetic cerrelations. See App jB 1| . 

In Fig. ^(a) mjj(p) is plotted for ^[u p ;x = 0] and *[ Zip j f j X = 0.1] for various choices of v. Finite size scaling in 
Fig.^(b) for x = 0.1 indicates vanishing long range order, mo — > 0, at p c = 3.3. It lowers the bound given previously 
by Ref. 0: at p c < 5. In Fig. ||(c) ??2q(£ _1 ) is plotted for $[u ex ; j; = 0] and ^[u ex ,v;x = 0.1]. Finite size scaling in 
Fig-l(d) indicates toq — > 0, at £ 1 = 0.3. In all the cases the GA (lines) works well. 
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Good agreement between GA and Monte Carlo is also seen in Fig.||(a) and (b), where S(r) is plotted for u p , and 
u ex respectively. Note how slow Sir) decays for = 0.3. By Fig.Q(b) mo(L) ~ L~ QA2 in this state. Exponentially 
decaying spin correlations are seen, both by Monte Carlo and GA, for u p with p > 3.7 and u ex with £ _1 > 0.4 Jl2| . 
Details of u sr , Eq.(|l6|), have strong effect on long range spin correlations. 

We use GA to extrapolate Monte Carlo calculations for S(r). In App.C 1 we find for exponential bond amplitude, 
u(r) — exp (—»"/£) and £ 3> 1, that S(r) decays exponentially with correlation length 



£ ex <x exp 



(23) 



For Gaussian bond amplitude, 
with correlation length 



u(r) = exp (— fir 2 ) with /i <C 1, we find in App.C2 that S(r) decays exponentially 



£s oc — exp (1 - x) — 



4/i 



(24) 



For Up, App.p 3] suggests vanishing long range order, too — ► 0, at p c < 3. 

Correlation lengths (|||) and ( p4| ) explain the slow decay of S(r) in Fig.||(b). It also indicate, that in the L = oo 
system, a small change in the ground state parameters brings an extremely sharp change in long range magnetic 
correlations. 



B. Superconducting order parameter 



The superconducting singlet order parameters are 



A i d = J2 Cs ' d ( f i) 



(25) 



where 



A« = flfUaibj - b iaj )/V2 



(26) 



The expressions of A's matrix elements are discussed in App.B3. By gauge invariance imposed by the Gutzwiller 
projector, (A s,d ) = 0. However, ty[u, v s (4y, x > 0] describes true s (d)-wave superconductors as seen by the singlet 

pair correlation function ((A")* ! Ag), a — s.d, in FigJ|. For v = v a , lim r ^ 00 ((A") t Ag ) ^ and ^ has (off-diagonal) 
long range order in A a . In contrast, the insulator states *&[u, Vi ns , x] and the "metallic" states ^[u,v metl x] have no 
long range superconducting order of either symmetry. 



V. EFFECTIVE HAMILTONIANS 



A. Magnetic energy and related parameters 



Magnetic order is driven by the diluted antiferromagnetic quantum Hcisenberg model 

H J = jJ2^r (27) 

(ij) 

Magnetic energy for x = 0: In Fig. Rfa) E mao (p), E mag (^) and E mag (Q) are plotted as a function of m^p), Tig(£) 
and to,q(Q) for u ex , u p and u g , Eqs. p3|), ( |14| ) and ([l^) respectively. In x = 0, all the three bond amplitude yield 
lowest magnetic energy of 

E a = -0.335 ± 0.0005 J /bond, (28) 

For Up, the optimal value of p is p op timai = 2.7, and m^(p = 2.7) = 0.105 ± 0.005. The ground state parameters of 
the Heisenberg model on an L = 40 lattice are : E(ground state)= 0.3347 J/bond and /^(ground state)=0.109 0. 
Table [| contains a summary of results for the optimal choice of parameters in all the classes. 



5 



Magnetic energy for x = 0.1: In Fig.|| E mag (p) and E mag (^) are plotted as a function of mg(p) and Wq^), for 
x = 0.1 and various choices of v from Eq. (|l7|). Within numerical errors, all states minimize H J at the same optimal 
parameters as for x — (Table ||). For u p , by Fig. [|(a) h yields local magnetization of m^(x = 0.1) = 0.08. For u ex , 
Fig. |](c) shows m 2 (x = 0.1) = 0.1. Thus we conclude that aside from the trivial kinematical constraints, the hole 
density and correlations have little effect on the magnetic energy at low doping. 

A better understanding of the properties of the optimal bond amplitude for H. J is gained by the average loop density 
defined below. From Eq. (^2j), a DLC contributes to tuq, Eq. (|2l|), it's number of pairs of sites, which share the same 
loop hence 



\AeA 7 / 

where l\ = ^2 ie x 1 is the loop length, and i 6 Aj. Thus L 2 TOq = .5(71", 7r) is proportional to the average loop length 
per site. 

The average radius of gyration of a loop is: 




(30) 



where 



w) j =rE('i-o'4i;("-'i) j . ( 3i ) 

i£A A ij6A 

with r^ m = X)ieA r «' anc ^ n A ^ s * ne number of loops in the DLC A. With Eqs. ( p9| ) and ( |30| ) we define the average 
density of a loop per site 

T= — J-?.. 32 

r 9 

The average loop density, T, is plotted in Fig.^(b), in the undoped case for all the bond amplitudes (fl3"|), 
and (p^[); and in the doped case for ^[tt p , w met ; a; = 0.1]. Comparison with Fig.^(a) shows that T is correlated 
with the magnetic energy. For vanishing mo, T converge to its value in the dimer RVB state, Eq. (|l2|), where 
r(dimer RVB) 9.6. This value of T is only slightly larger than T's value for an ensemble of DLCs, which include only 
configurations with two (or four) sites loops with dimer bonds. For such loops r* = 0.5 (or V2/2) and T = l\/r 2 g A = 8. 

The occurrence of loop lengths is interesting. In Fig ^9] we plot an histogram of the number of loops (fix), 
versus the number of sites on a loop (l\). The size of the lattice is L = 128, and u — iimf, which is derive from the 
Schwinger bosons mean field theory of TL J p5[. For all the bond amplitudes and lattice sizes we have checked n\(l) 
decays either algebraically or exponentially. 



B. Single hole hopping energy 

A single hole hopping in the antiferromagnetic background has been shown by semiclassical arguments Jl^Jl^] , to 
be effectively restricted at low energies to hopping between sites on the same sublattice: 

(ik)£A,B 

where i, k are removed by two adjacent lattice steps, and t' > 0. Unconstrained, the single hole ground state of H* 
has momentum on the edge of the magnetic Brillouin zone, in agreement with exact diagonalization of t — J clusters 
p9[ . Previous investigations have found that mier-sublattice hopping (the £-term in the t — J model) , is a high energy 
processes in the AFM correlated state [ ^8|Jl5| ] . We thus expect the same to hold even in RVB spin liquids with strong 
short range AFM correlations but no long range order. The primary effects at low doping may be to shift the ordering 
wavevector. 



G 



We denote by t' d (t' h ) the coefficients of second (third) nearest neighbor hopping terms. For t' h > t' d /2 the single 
hole bend minimum is at k m i„ = (±tt/2, ±vr/2), otherwise k m i„ = ±(0, ir), ±(ir, 0). Here we put t' h — l,t' d = 0.5. 

Results for the expectation value of H l are plotted in Fig.|l0|. The single holes hopping, Eq. (|33"|), prefers the 
metallic states v — v met over states with v — v s , va [ p^| - It also prefers longer range u(r) and thus actually enhances 
magnetic order at low doping. This is a type of a Nagaoka effect, where mobile holes separately polarize each of the 
sublattices ferromagnetically. 



C. The double hopping energy 

We consider Cooper pairs hopping terms 

H J ' =-j'lj2 A h A ^+ E 4 A ^"J ( 34 ) 

\ijk ) 

Calculation of Tt J matrix elements is discussed in App.|B 3|. The first term in TL J is derived from the large U Hubbard 



Term (a) is a rotation of the singlet pair. It 
is a parallel translation of singlets. It prefers 



model to order J' = t 2 /U jl5|]. It includes terms (a) and (b) in Fig. 11 
is positive for v s fT^] and hence prefers Vd over v s . Term (c) in Fig. 11 
superconductivity with v = Vd or v = v s over metallic states with v — v met ppf . For x = 0.1, Ti. J is minimize by Vd- 

In Fig. ^ the ground state energy E p h of (Q) is plotted for v = Vd, u — u p and u — u ex . a;=0.1, and the size of the 
lattice is L=40. The variational energy is minimize at p = 3.35 and = 0.35, for u ex and u p , respectively. In both 
cases, by the finite size scaling of FigJJ(b) and (d), it indicates vanishing tuq at L — > oo. Thus, Cooper pair hopping 
drives the groundstate toward a spin liquid phase! 

The Gutzwiller approximation fails to predict this effect. According to the GA, the minimum of the double hopping 
energy roughly coincides with the minimum of the magnetic energy ((Ti. J )). This is understood by (see App. ^): 

{S+Sj) GA = -{a\b])l A , 

where i G A, j G B. The GA agrees with Monte Carlo results for matrix elements of long range pair hopping. 

The matrix element of (d) in Fig.[TT|, and (njuj) also drives the groundstate toward a spin liquid, and prefer 
superconducting over metallic states JL2[. These terms are excluded due to relatively large thermal noise. 



VI. A RELATION BETWEEN SUPERCONDUCTING T c AND LOCAL MAGNETIZATION 

Since Ti. J is the effective model which drives superconductivity it produces phase stiffness, which in the continuum 
approximation is given by 



H J ' 



^ J d 2 x(V0 4 ) 2 (35) 



The stiffness constant Vq can be determined variationally from the doped RVB states. Imposing a uniform gauge field 
twist on A, Aij — > Aijexp(i(xi + Xj)<j)/2L), (H J ) becomes, to second order in <p/L, 

E P h- — 

V = ^ = 2J' (<A^A , £ ) + (Aj . A ,_ fi ) + <A^A £+ ^>) (36) 
Following Ref. ||, at low doping for the square lattice Vq is roughly equal to T c 



In Fig.(|l3]) we show our main result: The staggered magnetization mo(p) for TL J + H J is plotted against the 
superconducting to magnetic stiffness ratio Vo(p)/J for different doping concentrations x = 0.05,0.1,0.15, v = Vd, 
and u = u p . The actual free parameter in the graph is J'/ J, from which mo and Vq are determined variationally. 
Two primary observations are made: (i) The local magnetization is sharply reduced at relatively low superconducting 
stiffness (and T c /J). (ii) The relation between mo and Vq/J appears to be independent of x. 

For u ex , Eq. @, it requires Vq{£)/J = 0.49 for (H J + H J ') to be minimized at ^ x = 0.3. By FigJ|(b) this leads 
to m^ = °° = 0. 
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VII. SUMMARY AND DISCUSSION 



In this paper we used extensive Monte Carlo calculations to study properties of hole doped RVB states. We 
found that an effective model which include Heisenberg and pair hopping terms is consistent with the experimental 
connection between superconductivity and reduction of local magnetic moment. Within checked variational options 
we showed that the properties of the model are independent of particular choice of parameters for the state. Gutzwiller 
mean field approximation for magnetic correlations was found to agree with Monte Carlo calculation, and used for 
analytical extrapolation of numerical results. We showed that long range magnetic correlations in RVB states are 
extremely sensitive to variational parameters. We found that the average loop density is well correlated with the 
magnetic energy. We conclude this paper in several arguments and insights regarding our results. 

Magnetic energy and long range magnetic correlation: Note the contrast between correlation lengths ( ^3|) and (|2~i|), 
and the "shallowness" of the minima of the magnetic energy in FigJ^. It imply that a very weak pair hopping term 
in the Hamiltonian causes a dramatic change in long range magnetic correlations. 

Magnetic energy and loop density: A comparison between loops (a) and (b) in Fig|l4| shows that large amplitude 
(Qa) of DLCs with "denser" loops enhance the probability to find nearest neighbor sites on the same loop and reduce 
the magnetic energy. 

The loop density shows that the optimal bond amplitude is determined by an intricate balance between mo and r g . 
This relates quantum spin fluctuations to the average loop density of the ensemble. 

Effective model for doped system: Ti. J + Ti. J describes the low energy physics of the lightly doped Cuprates. As 
the lattice is doped, its variational ground state is a d-wave superconductor, with a sharply reduced local magnetic 
moment. The model includes built-in pairing. Such a model is supported by the existence of a pseudogap in the 
normal state of the high Tc materials. 

Relation between phase stiffness and local magnetization: Because of finite size uncertainty, TOq = °° in Fig.jfilf) is an 
upper bound on the thermodynamic local magnetization. A sharper reduction of the local magnetization occurs if: 

(a) The GA result of App. |C 3| , mo(u r — r~ 3 ) = 0, is correct to the discrete lattice. In that case mo vanishs already 
at V /J > 0.2. 

(b) In finite doping the optimal bond amplitude for Ti. J 



TL J decays exponentially. In that case mo vanishes for 
Vq/J > 0.5. Variationally, we can not rule out this possibility. 

In both of these cases there is a qualitative agreement with the doping dependence of the local magnetization and T c , 
as measured by Refs. HQ. 

Useful conversations with C. Henley, S. Kivelson and S-C. Zhang, are gratefully acknowledged. MH thanks Taub 
computing center for support. A A is supported by the Israel Science Foundation and the Fund for Promotion of 
Research at Technion. 



APPENDIX A: THE FERMION PART OF THE DOPPED RVB STATE 



The fermion part of \^f[u, v, x]) is 

\9(x)) f =P G (N h )exp 
where Nh = xL 2 . We write this state as 



E v«ftfl 

i£A,j£B 



E ^flf] 



|0) 



(Al) 



N h /2 

i*(*)>/=E c (7)n/Uj°> 



fc=i 



(A2) 



where 7 is a distinct configurations of holes sites: 

7 = {(ik G A, jk G B)} k 2 =1 : Vfc < k ik < %h> jk < jk> 



From Eq.flAfl) 



l*(W/ = il>ai/Ji 



2 ■ 7 
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where ai S {ik}^={ 2 , Pi <= {jfc}^i 2 - For each 7 we commute pairs of operators, without any affect of sign, and order 
A holes operators in an increasing order of their site index 

y CT 2 2 2 " l 2 J 

where a — a (Nh/2 — > Nh/2). Commuting the _B operators we get 

C(7) = detF (7) 

where V is an x Nh/2 matrix with 

When 1; connects same-sublattice sites the determinant is replaced by a Pffian Q . 



APPENDIX B: EXPRESSIONS FOR EXPECTATION VALUES 



1. Alternative calculation of magnetic correlations 

We use the operator identity Eaj 



At. A,, 

*7 y 



S^Sj 



1 



(1 - /*/«)(! - /)/;), 



where Ajj = // n(aibj — bidj)/\/2, to express (SjSj). 
show that for i E A, j E B and u(r) > V?' 



SiSj (7,A 7 ) 2 



1 x 

4+2 



"In Its i 



(i,n)ec 



(i, j) is a bond in a (A 7 ) 



i, j <E 7 
otherwise. 



a (A 7 ) is the set of forward bonds in A 7 , of the sites on sub-lattice A. 

We demonstrate Eq. (p3) for half filled lattice ({fjfi) = 0). With M tJ = (a»6j - Mi), 



M <i M+|0)=2|0) 



(Bl) 



(B2) 



(B3) 



and hence 



A4My|*[tt]) = ^ 



J 



n uin i a y/ 



. ((,n)ec«j 



(B4) 



where |o!y) is a valence bond state, with (i, j) G a^ . 

The term in the square brackets requires further explanation. From Eq. (B4), for any pair (l,n) E onj : / =/= i, 

\ctij) = Mj.Mij\/3), where (i, j), (l,n) $ (3, (i,n) E (3, and otherwise (3 = a , see Fig.(^5|). In each \0) 



carries a factor uijUi n . Eq. ( JB3| ) indicates an additional option to get o^, from onj. 

Taking the overlap with we get the matrix element which is expressed in terms of Eq.(B2). aiy represents 

the Ket. A possible definition of the bonds of the Ket is the forward bonds of the sites on sublattice A. 
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2. Matrix element of single hole hopping term 

Fig.|l6| describes the effect of a single hole hoping term on a hole-pair configuration. Using definition (p"l|), for 

i, k G A 



dctV(7 fc ) uu 



fJfk(ala i + blb i )( 7 ,A 7 )=l ^vhYu kl 



^-s kl if i £ 7, fc ft 7 







otherwise 



(B5) 



where i g" 7^., fc G 74., and otherwise 7^ = 7; (fc, Z) € A 7 is the forward bond of k (i.e. originated in the Ket); and 
Ski = ±1 comes from reordering the fermion operators. Relation (B5) is simplified using p2[ 



detV(7fc) s ki = detV (7, i -> k) 



(B6) 



with 



V (7, 



v(k,j p ) i r = i 
V(l)rp otherwise 



(B7) 



3. Matrix elements of the double hopping terms 



For u r > Vr we express (A t, A. ,,-), where A is given in Eq. (26). Aj. ; creates a singlet bond. Ay creates a pair of 
holes. With the results of App.B 1, for i E A, j E B 



( 



s dcty( 7a ) 
2u kl detV(7) 







E 



2ui 



. (r,n)6a(7) 



if (k,l) G a( 7 ), 
i E 7 if i 7^ k, 
j G 7 if j 7^ Z. 



otherwise 



(B8) 



Where s = —1 if i = k exclusive-or j = I and 1 otherwise, V(-y a ) = V(^,i — > k,j — > Z) is defined like Eq. (|B7j), with a 
possible replacement of a row and a column, and 0(7) is the set of forward bonds of A sub-lattice sites in A(7). 



APPENDIX C: THE GUTZWILLER APPROXIMATION. 



The Gutzwiller Approximation amounts to dropping the projector T'(x) in definition (^) and setting |\&[it, u; 2;]) 
\ip[yu, zv\) ~ \yu) ® \zv). The constants y = y(u) and z — z(v) are determined by global constraint equations 



{a\ ai ) = (bth) = (1 - x)/2, 
(flfi) = x 



(CI) 
(C2) 



for y, z respectively. In this section (...) = (^]...|^)/(^1^). 

\yu) is a Schwinger bosons mean field wave function Jl5| , on which we preform the Marshall transformation [ p"5| 
aj — > — bj, bj — > a,j, j G B. Hence 



\yu) - exp y^iMa] + b\b]) |0). 



(C3) 



Operators are transformed accordingly, for example Sj — > — Ojbj for j G B. 
From Eqs. del) and @ 



1 ^1 • xf 



(C4) 
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Whereas (¥|S+Sr(l - / l t /0 2 |*)/(*l*> = (1 - x)/2. Thus we use 

(l-x/3)- l {StSr) (C5) 

as the GA for the long range magnetic correlations and mo in the doped RVB state. Empirically we omit the 
(1 — x/3) -1 factor in the estimates of magnetic energy. 
Using the extended Wick theorem |H| , for i e A 

(q + 9 -v = / -<«l & i«] 6 i) = -(«!«]) (^^') = -4- J e B (CM 
[ * i> \ (athbfa) = (a\a 3 )(b)h) + Sy/2 = a* + %/2 j € A K ) 



where we used, for example, (a\bj) = (aj&j) = 0. 

Expanding «k = Ylj e lkj Uoj, Pk = X), e ^ Poj an d a similar expression for Ok pl| 



and 



The constraint equation ( |C1| ) becomes 



We consider three cases: 



1 - y^u£ 



= T^h- (08) 
1 - 



La2i-it£ 2 ■ (C9) 



1. Exponential bond amplitude 

We calculate spin-spin correlation function for 

u(r) = 

with £ >> 1. 



exp (— r/£) r bipartite 
otherwise 



= J J d 2 r exp (-r/£) exp (-ik • r) = 2tt£_- 1 / (k 2 + £ 



-2\ I 



where MBZ= magnetic Brillouin zone. Eq.(C9h becomes 



k e MBZ 

otherwise 



2y 2 /• d 2 fc _ Try 2 

~F J mbz (fc 2 + C- 2 ) 3 - (2^- 1 y) 2 ~ IF" fc 3 - (2^-^) 



1 - x 



(CIO) 



(OH) 



(C12) 



where we multiplied the left side in 2 to account for the integration over the complete Brillouin Zone. In all our 
calculations we took the continuum limit (lattice constant— >0), where the upper bound of the integration — > oo. This 
appro xima tion works very well for slow decaying bond amplitude jl2| . 
Eq.(|Cl~2|) gives || 



1 



3(2tt)^ U 



In 



Try 



-^2V3arctan 
'2ttV 



c8 



23 



1 

2 In 



Hi 



7T 2 y 2 ^ 4 
Try 



i 



(C13) 
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The argument of the last logarithm has to be sufficiently close to zero for Eq.(C13) to be satisfied. Therefore 

1 



2tt£ 2 



Hence we can neglect in the left side of Eq.(C13) all terms but the last log. Consequently 



1 



1 3 ^ 2 

2*4 1 l ~ — eX P 



2^ 



(l-x) 



3< 2 



Eq.(C6) becomes 



a T = 2 / d 2 ktr k exp (ik • r) = — -r / 

7mbz 2 n j 



- {2-Kfe 
kJo (kr) 



V 



(k 2 + £- 2 f - £- 6 V 



(C14) 



(C15) 



(C16) 



where Jo is the Bessel function p2j. Since the integrand in Eq.( C16| ) vanishes as k — > 0, we can replace Jo with its 
approximation for fcr 3> 0. Expanding the denominator to first order in k 2 



V 



dk-. 



Vk 



Gwy/rTr^ 2 J k 2 + a 2 
where a 2 = (1 - X>)/(3£ 2 ). In the definition 



[cos (fcr) + sin (kr)] 



V 



;Y 



(C17) 



' /s 2 + a 2 



Yo = i?e Yi + JmYi . Let us consider the Integral 



dz- 



■ exp (izr) 



exp (ikr) 



dk , y~ exp (i/cr) 
fc 2 + a 2 y ' 



where the close contour encircles the upper half of the complex plane. The part of the contour along the negative real 
axis is 

Vk f°° i/fc 7 

dkj 2 o ex P {ikr) —if dk' '— % ex P (—ik'r) = iY^ 

where we substituted k' = e~ m k. Hence Y2 = Y\ +iY* = (l + i)(ReYi +ImYi) = (l+i)Yo. Using the residue method 
for Y 2 : 



exp 



-ry/(l V)/(3e) 



and with Eq. (C15) we find for the correlation length of the spin-spin correlation function, £ e 

/ ^ 1 
~ utt; — = 23 ex P 



4(1-1?) 



3tt ^ 



(C18) 



(C19) 



2. Gaussian bond amplitude 

We calculate spin-spin correlation function for 

exp (— fir 2 ) r bipartite 







otherwise 



For k e MBZ 
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TT / k 2 \ 

Uk= ^ CXP { J 



Eq.Q is 



with the solution 



1 -a; 



2 A*' 



-(1-z) 



2/i 



7T^ 



Calculation of a r is identical to the exponential case. Substituting in Eq. ( C17 ) a 2 = 2p(l — Q) 

exp (- r j2n{l-G)\ 
<J T (X -= 

and hence the spin-spin correlation function decays exponentially with correlation length 

1 1 



y/8fJL(l-g) 



exp 



(1-x) 



4/i 



(C20) 



(C21) 



(C22) 



(C23) 



(C24) 



3. Power Law bond amplitude 



For the bond amplitude 



I — T r bipartite 

otherwise 



(C25) 



we show, in the continuum limit, that for < e < eo, S(tt, 7t) is finite and hence mo = 0. 

Calculations of the GA on lattices of size L < 512, show that for any e, the spin-spin correlation function calculated 
with function (C25), decayes slower than with u — 1/r 3 . This suggests that toq = for u — 1/r 3 . 



s(n, tt) = £ i (s+ 5 -) i = £ < + £ + 1 = -i y: « +pd+1 



where we used Eq.(OJ). For k e MBZ p 



lik = 27re e 



2 -efc 



(C26) 



From Eqs.(|C7|) and (C8) e efc , and <7k *— e 2efc . Hence 5(7r, tt) might diverges only if there is ko, such that 



A; — >oo 



, 1 — ae 



-2ak\ 



' k—kn 



0, where a = a(e) = (27rye 2 ) . Therefore if a < 1, 5(7r, tt) is finite. 



Eq.(p9h for y is 



which becomes 



a f°° „ fce" 2efc 
7r Jo 1 — ae zeK 



a r ~dk^ 



1-x 



(C27) 



(C28) 



P =i 



The right side of Eq.(C28 ) is 3/ independent, and increases with e, hence a increases with e. Therefore a(eo) = 1. For 
a = 1, the left side of Eq.(]C2§) is tt 2 /6 and 
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and for e < eg, S(n, n) is finite and hence mo = 0. 



state 




Eo(x = 0) [J/bond] 


Eo parameters 


mf)(x — 0) 


u ex , Eq. ( 


lij) 


-0.335 ± 0.0005 


C 1 = 0.17 ±0.005 


0.125 ±0.005 


u g , Eq. ( 


15 


) 


-0.335 ± 0.0005 


Q = 0.014 ±0.001 


0.12 ±0.005 


Up, Eq. ( 


13 


) 


-0.335 ± 0.0005 


p = 2.7 ±0.05 


0.105 ±0.005 


UMF 


-0.3344 ± 0.0002 




0.087 ± 0.005 



TABLE I. Minimal magnetic energy (Eo), optimal choice of parameters, and square of the ground state magnetization (mo), 
for various bond amplitudes in the undoped, x = 0, case. The size of the Lattice is L = 40. umf is derived from the Schwinger 
bosons mean field theory of T( J . The ground state parameters of the Heisenberg model on an L = 40 lattice are : E(ground 
state)= 0.3347 J/bond and mo(ground state)=0.109 (Tr). 
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FIG. 1. A bond configuration in the doped RVB states ^[m, «]. Solid (empty) circles represent spins (holes) with bond 
correlations Uij (vm). 
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FIG. 2. Exact and numerical Monte Carlo results of magnetic energy and correlation length of the spin-spin correlation 
function, vs. hole doping- a; H|. Dimer v and u, Eq. (|l2|), is used on a 2 x 64 ladder. 
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FIG. 3. v = v^ et . £ is centered at Fermi pockets around k m i„ = (±7r/2, ±7t/2). Within the pockets u me t = ±1, otherwise 
Vmet = 0. Note that v me t has d x 2_ y 2 symmetry. 
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FIG. 4. (a) The local magnetization squared, m,^(p), of doped and undoped RVB wavefunctions, versus the variational power 
p, defined by Eq. ( |l3| ) . Lattice size is 40 x 40, and in the doped case hole concentration is 10%. Results agree with the Gutzwiller 
approximation (lines). The hole bond parameters v are defined in Eq. (|l7|) . is weakly dependent on v. (b)Finite size scaling 
of mo(L) for p = 3.3 which indicates vanishing local magnetization at L — > oo. (c) mo(£) versus the variational correlation 
length £, defined by (0). (d) Finite size scaling of mo(L) for £ 1 = 0.3 which indicates vanishing local magnetization at 
L — > oo. 
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FIG. 5. (a) Calculations of spin-spin correlation function (S(r)) with Monte Carlo and GA (lines) for undoped states with 
Eq. The size of the lattice is L — 128. (b) Calculations of S(r) using u — u ex , Eq. (y). Note how slow S{r) 

decays. In both cases there is a good agreement between Monte Carlo and GA, and S(r)) weakly effected by doping. 
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FIG. 6. The singlet pair correlation function. ((A" Aq). A" defined in Eq. J25|), and a — s,d. In this graph u = 1/r 3 ' 3 . 
ty[u, v a \x] has superconducting (Off-Diagonal) long range order only of symmetry a. ^[u,v m et',x] has no ODLRO in either 
symmetry. 
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u g ,d=0 
u ,5=0 

ex ' 

u ,5=0 

p ' 

u n ,h=0.1, v 



A 



met 



(a) 



5lb • a A 



■a 



- d 



\ © 

B 



0.1 



0.2 



m, 



FIG. 7. (a) Magnetic energy (Em afl ) versus local magnetization squared mj. Variational parameters are : for x = 0(left y 
scale), u p (p) Eq. (p"3|), u ex (£) Eq. (^4[), and u g (Q) Eq. (|l5|). For a; = 0.1 (right y scale) v me t and u p . lattice size is L = 40. 
The optimal parameters for each u appear in Table [IJ The minimal magnetic energy for x — is —0.335 ± 0.0005J/ bond, (b) 
The average density of a loop per site V, Eq. (^), versus nig for variational cases as in (a). In all the cases F is correlated with 

Emag . 
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FIG. 8. Magnetic energy (Emag), for u v Eq. ([H^), and u ex Eq. ([wj), versus local magnetization squared mo, using various 



hole distributions from Eq. 



The density of holes is x = 0.1 and lattice size is L = 40. E mag is weakly dependent on 



inter-hole correlations. For u v (u ex ), E mag is minimized at mo(p) ~ 0.08 (0.1), p = 2.7 (£ 1 = 0.17). 
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FIG. 9. The average number of loops (n\), versus the number of sites on a loop (lx). The state is with u — umf, which is 
derive from the Schwinger bosons mean field theory of H J . The size of the lattice is L = 128. For 1x^,130, ni oc l x 2 ' 7 . 
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FIG. 10. The single hole hoping energy Ej lov , Eq. (P3|), versus local magnetization squared m\ for metalic hole distributions, 
and spin bond amplitudes u p and u ex . The density of holes is x = 0.1 and lattice size is L — 40. In H l , t' h = 1, t' d = 0.5, 
and the single hole band's minimum is at k m i n = (±7r/2, ±7r/2)). The single hole hopping prefers longer range u(r) and hence 
higher local magnetic moment. It also pefer metallic states over v — v s ,Vd [fl^l - 



T 

l> ■ 



o - -o 



(a) / (b) (c) (d) 

FIG. 11. (a), (b), and (c) are the terms of H J , Eq. (|34j). Dashed line and empty circles = A. Term (a) is a rotation of a 
singlet pair, it distinguishes between s to d wave superconducting order parameters. Term (c) prefers vm s ) over metallic states 
with Vmet- Term (d) dependency on the variational parameters is similar to that of (c), it is excluded due to thermal noise. 
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FIG. 12. The expectation value of H J , Eq. (^), E p h, versus mo, for u v and u ex . In contrast to the magnetic energy Fig. 
E p h prefers a vanishing mo at L — > oo. Note how similar the graphs are for u ex and u p . 
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FIG. 13. The relation between thermodynamic local magnetization mj 00 and superconducting phase stiffness Vb (related 
to T c , see text), u = it p , Eq. (|l3|). J is the Heisenberg exchange energy. The points are considered upper bounds on mo, which, 



for Up, may even vanish for Vo/J > 0.2. For u ex , Eq. (0), mo vanishes for Vo/J > 0.5. 
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dilute 



(a) (b) 

FIG. 14. Two kind of loops: (a) "Dense loops", which fully cover small regions of the lattice, and many nearest neighbor 
pairs. Bond amplitudes u which maximize the weight (Ha) of loop configurations with such loops minimize the magnetic energy, 
(b) "Dilute loops", which contributes very few nearest neighbor bonds to the magnetic energy, Eq. (p^). Loop (a) is denser, in 
the sense that it covers more sites on roughly the same "area"= (r^) 2 , Eq. (|3l|). 
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i i;i _tk)i 

1 o • — • 1 o • — « 



cx(Ya) «(Y) 

FIG. 16. The operator ///fe(a^ai + turn a hole pair configuration 7*,, 0(7*,) (left), with k g 7*,, 6 a(7fc), to the 
right configuration 7, 0(7) with 167 and (fe, ?) 6 0(7). In f} fk{a\.ai + &J.&j)|\E'(, this configuration has the coefficients ua, and 
detV(7 fc ). 
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